Lipid membranes with an edge 
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Consider a lipid membrane with a free exposed edge. The energy describing this membrane is 
quadratic in the extrinsic curvature of its geometry; that describing the edge is proportional to 
its length. In this note we determine the boundary conditions satisfied by the equilibria of the 
membrane on this edge, exploiting variational principles. The derivation is free of any assumptions 
on the symmetry of the membrane geometry. With respect to earlier work for axially symmetric 
configurations, we discover the existence of an additional boundary condition which is identically 
satisfied in that limit. By considering the balance of the forces operating at the edge, we provide a 
physical interpretation for the boundary conditions. We end with a discussion of the effect of the 
addition of a Gaussian rigidity term for the membrane. 

PACS numbers: 87.16.Dg, 46.70.Hg 



I. INTRODUCTION 

A closed lipid membrane is described remarkably well 
by a geometrical liamiltonian. This hamiltonian is con- 
structed as a sum of the scalars, truncated at an appropri- 
ate order, which characterize those features of the mem- 
brane geometry which are relevant. A term quadratic in 
the extrinsic curvature provides a measure of the energy 
penalty associated with bending jj], ||, ||, |I| |); any in- 
trinsic tendency to bend one way and not the other is 
captured by a term linear in the extrinsic curvature [[3) . 

The shape equation determining the equilibria of this 
membrane is a fourth-order non-linear elliptic PDE of 
the form V 2 K+ 'A' 3 ' = 0, where V 2 is the Laplacian on 
the membrane, K is the sum of the principal curvatures, 
and by 'A 3 ' we mean a cubic polynomial in these curva- 
tures Here, we would like to examine the boundary 
conditions which must supplement this equation when 
the membrane possesses a free edge. The energy cost 
associated with this edge is, to a good approximation, 
proportional to the exposed length. During the forma- 
tion process, material will either be added to the edge or 
the edge will heal itself so as to form closed structures. 
However, unlike a simple soap film with a free edge, a 
lipid membrane does possess stable (cup shaped) equi- 
libria Q . An examination of such structures is a small 
but important step towards a fuller understanding of self 
assembly. 

Our primary focus will be on the boundary geometry. 
We have a surface with a boundary and a certain energy 
penalty associated with it, a well defined problem in clas- 
sical field theory. The boundary conditions are identified 
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by demanding that the energy be stationary for arbitrary 
deformations of the edge geometry. In distinction to ear- 
lier analytic treatments of this problem, we will relax the 
assumption that the membrane geometry be axially sym- 
metric Q . This is important not only for conceptual rea- 
sons. Generally, there will be no priveleged parametriza- 
tion such as that tailored to axial symmetry; in an axi- 
ally symmetric geometry the edge itself is simply a circle. 
We find that there are three boundary conditions. As we 
will demonstrate one of these conditions, involving three 
derivatives of the embedding function, is satisfied identi- 
cally in the axially symmetric limit. Therefore this limit 
cannot be considered as a reliable guide to the general 
case. 

While the variational approach does capture the ge- 
ometrical nature of the boundary conditions, the phys- 
ical interpretation of these conditions still needs to be 
clarified. Ideally, one would like to interpret them in 
terms of the balance of forces operating at the edge. To 
do this in a way which does justice to the geometry, we 
identify the conserved Noether currents associated with 
the intrinsic translational invariance of the configuration 
|j. The three, apparently unrelated, boundary condi- 
tions are now cast in terms of the three components of a 
single vector identity on the edge. 

We finish with a discussion of the effect of a Gaussian 
rigidity term on a lipid membrane with edges. Whereas 
such a topological term does not alter the bulk shape 
equation, we show that it does modify the boundary 
conditions which apply to it in a way which will have 
consequences in the bulk. This extension is relevant in 
topology changing processes. 

The outline of the paper is as follows. In Sect. II we 
consider the simple example of a surface tension domi- 
nated membrane. This allows us to establish our nota- 
tion and to derive the boundary conditions in a simple 
context. In Sect. Ill, we derive the boundary conditions 
at the edge for a lipid membrane. We then specialize to 
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axially symmetric configurations to compare our results 
in this limit with previous work on the subject. In Sect. 
IV, we consider the balance of the forces operating at the 
edge and we show how they are related to the bound- 
ary conditions. The effect of a Gaussian rigidity term to 
the membrane energy is the subject of Sect. V, where 
we obtain the appropriate modifications in the boundary 
conditions. We end with some remarks in Sect. VI. 



II. SURFACE TENSION PITTED AGAINST 
EDGE TENSION 

It is useful to examine first the simpler situation in 
which the membrane physics is dominated by surface 
tension, such as a soap film with a free edge. Let the 
membrane surface E have an area A, with boundary C 
of length L, and the tension in the membrane bulk be a 
constant /i, and that on the edge a. The energy is then 
given as a sum of two terms 

F = liA + <jL. (I) 

Surface tension tends to decrease the membrane area; 
line tension to decrease the length of the free boundary. 
Without some further refinement this model does not ad- 
mit stable equilibria. Suppose a hole is punctured in the 
film, then depending on its radius, either the hole will 
close healing the film, or grow and destroy it. An unsta- 
ble equilibrium clearly exist when the radius is tuned to 
coincide with a critical value r c . On dimensional ground, 
one world expect r c rj fi/u. (Wc will ignore stability here 
as our interest in this model is only as a point of reference 
for a lipid membrane.) 

The membrane surface E is described by the embed- 
ding X in three dimensional space R 3 as x = X(£ a ), 
where x are coordinates for R 3 , and £ a coordinates for 
the surface (a, b, ■ ■ ■ — 1,2). Its edge C is embedded in 
turn as a curve on E as £ a = Y a (s) which we parametrize 
by its arclength s. We can now cast F as 

F = M / d 2 tV7 + vfds. (2) 

Here, the metric induced on E is given by 7 ao = e a • e^, 
where e a := <9 a X are tangent to the surface, 7 = det7 a b, 
and dA — ^7 g? 2 £. Note that we can also consider the di- 
rect embedding of the edge C in i? 3 , via x = Y(s), where 
Y = X(V a (s)). The tangent to C in R 3 is equivalently 
expressed in either of two ways: t = e a t a , where t a = Y a ; 
or t = Y, where a dot denotes a derivative with respect 
to arclength s. 

The energy is a functional of the embedding X of E 
in R 3 . There is no need to vary the edge embedding Y a 
independently: the Y a are fixed by the constraint that 
the two embeddings agree, Y = X, on C. 

Equilibrium configurations arc those at which the en- 
ergy Eq.(|l|) is stationary. To derive the equations de- 
scribing the equilibrium configurations of a membrane 



described by this model, we first consider a variation of 
the embedding X of the membrane X — > X + 5X. We let 
n denote the unit normal to the surface E. We decom- 
pose the displacement with respect to the spatial basis 
adapted to E, {e a ,n}, as, 

SX = $ a e a + $ n . (3) 

We now have that the induced metric varies according to 
(see e.g. §) 

Sxjab = 2K ab $ + V„$ 6 + V b $ Q , (4) 

where K ab denotes the extrinsic curvature tensor, 

K ab = e b ■ d a n , (5) 

and V a is the covariant derivative on E compatible with 
"fab- The derivative terms in the variation of ^ ab are its 
Lie derivative along the tangential vector field, $°. The 
variation of A is 

S X A = / dA [K $ + V a $ a ] 

= J dAK $ + <j> dsl a $ a . (6) 

The mean extrinsic curvature, K = K ab -f ab . The second 
line follows from the preceding one using Stoke's theo- 
rem. Here l a is the outward pointing normal to C on 
E. Only the normal projection $ of the variation plays 
a role in determining the bulk equilibrium of the mem- 
brane. This is true generally, regardless of the model. 
In this particular model, however, there is no boundary 
term associated with the bulk normal displacement <&. As 
we will see, this is not generally true — a happy accident 
when the energy is truncated at the area term. On the 
other hand, the tangential bulk variation $ a always gives 
only a boundary term. This is a consequence of the fact 
that a tangential deformation corresponds in the bulk to 
an infinitesimal reparametrization of the surface. There 
is however a physical displacement of the boundary. In 
fact, the boundary contribution to Eq.(|^) is easily identi- 
fied as the change in the surface area of E under a normal 
deformation of its boundary, SY a — (l b & b )l a , which at 
each point is directed along the tangent plane of E at 
that point. The projection of <& b onto the edge C itself, 
t a Q a , does not contribute. 

Let us turn now to the variation of the edge embed- 
ding Y induced by the bulk variation <5X. It can be 
decomposed with respect to a basis adapted to both em- 
beddings, X and Y, given by {t,l, n}, where 1 = e a l a . 
Thus at the edge we set 

5Y = <j> t + V 1 + $ n , (7) 

where the edge and bulk components are identified by 
continuity, ip = l a $ a and <j> = t a & a . Modulo a diver- 
gence associated with a reparametrization of the bound- 
ary which involves the tangential component 4> that we 
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can safely discard, we have for the variation of the in- 
finitesimal arclength 



5 Y ds = ds + K\\ 
where we have used the fact that 

t = k! — Khxi. 



(8) 



Here n is the geodesic curvature of C associated with its 
embedding in S, and we have defined Ku = K ab t a t b . 
The unconventional minus sign in the first term of Eq. 
(^) comes about because 1 is the outward normal to C on 
S, i.e. i a = -k l a . 

The corresponding deformation in L is then given by 



5 Y L = (h ds (k V + K\\ • 



(10) 



Summing the two contributions (|^) and ( |ic| ) to the vari- 
ation of the energy F, as given by ([!]), we find 

S X F = fx / dAK $ + f ds [(fx + a K)tp + a K\\ <fr] . 

(11) 

The bulk equilibrium is a minimal surface unaffected by 
the boundary, satisfying K = 0. On the boundary, the 
projections along the normals to the edge, tp and rep- 
resent independent deformations, so that stationarity of 
F requires the vanishing of the corresponding coefficents. 
We thus read off the two boundary conditions: 



a k + fx 
o-Ku 



(12) 
(13) 



The first tells us that the geodesic curvature of the edge 
as embedded in the membrane is constant. The second 
simply enforces the vanishing of K\\ at the edge. Note 
that the completeness of the basis {t, 1} of tangent vec- 



tors on £ at C, 



r 



t a t b 



l a l , permits us to ex- 



Note that had we N sheets conjoined on a single edge, 
Eq.(|l4|) gets modified in an obvious way: 



N 



(15) 



where 1^ is the vector normal to the edge which is tangent 
(9) to the i th sheet. Eq. (|l5|) provides a generalization of 



the Neumann rule for soap bubble clusters at a Plateau 
border [[to) to accomodate line tension on the edge. A 
simple application is considered in pi). 



III. LIPID MEMBRANE WITH AN EDGE 

A lipid membrane is modeled by a phenomenological 
energy quadratic in the extrinsic curvature of the surface. 
Let us write this as 



F b = J dA F{l a \K ab ) 



(16) 



i.e. T depends at most on the extrinsic curvature, and 
not, for example, on its derivative V a K bc . In particu- 
lar, we will focus on the model described by the Hclfrich 
energy density 



T = a(K - K ) 



(17) 



The spontaneous curvature A'o is a constant, as is the 
bending rigidity a. The constant [i is interpreted here 
as the Lagrange multiplier implementing the constraint 
on the membrane area. We will discuss the addition of a 
Gaussian rigidity term in the next section. 
The energy of the bulk and the edge is 



F = F b + aL. 



(18) 



The shape equation describing the equilibrium in the 
bulk, which is derived from the extremization of the en- 
ergy © 



press the mean curvature at the edge as A = A,| + K±, a[-2V A + 2(A - K )K + (K Q - A ) A] + \x A 
where K± = K ab l a l b . Thus modulo the bulk equilibrium 
A = 0, the boundary condition (|l^) can be alternatively 
expressed as K± = 0. The only potentially non van- 
ishing component of K ab on the edge is the off-diagonal 
component, K\\±_ — t a l b K ab . 

For this particular model our approach has been heavy- 
handed; the boundary conditions we have written down 
are an elaborate way to express the simple vector identity 



0, 
(19) 

is well known |7J . The structure of this equation has been 
discussed in detail elsewhere Q , where a novel derivation 
is also provided. The scalar curvature 1Z appearing in 
Eq.(|l9|) is related to the extrinsic curvature through the 
Gauss-Codazzi equation, 



11 = K 2 — K ab K ab 



(20) 



at = /xl , 



(14) 



which equates the change in the tension over the interval 
As along the edge, crAt, to the force due to surface ten- 
sion acting on the edge, /xlAs. The apparent mismatch 
in counting (three versus two) is accounted for by noting 
that the projection of Eq.(|l4|) along t is an identity. For 
higher order models, as we will see, this projection will 
not be vacuous. 



Under a tangential deformation of the surface, <5mX = 
$ a e a , the energy density transforms as a divergence 
which is transferred to the boundary, 



S l{ F b = j> dsTl a <b a . 



(21) 



This is because the local scalar energy density T trans- 
forms as 



(22) 
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The details of T do not enter. Note that Eq.(|21|) agrees 
with the corresponding expression for the area with T = 
1. As before, this boundary term induces a source into 
the boundary Euler-Lagrange equation. For an edge with 
a line tension a, we get the first boundary condition, due 
to a deformation along the normal 1, ij), 



an + T = , 



(23) 



where we have used Eqs. (|10|) and (|2lj). This should be 
compared with Eq. (|l2|) to which it reduces if T — /i, a 
constant. This boundary condition relates the geometry 
of the edge to the extrinsic curvature of the membrane 
evaluated at the boundary. 

We now examine a normal deformation of the surface 
E, S±X. = <f> n. The shape equation ( |l9|) determining 
the local membrane equilibrium is obtained by demand- 
ing that the energy be stationary with respect to normal 
deformations of £ which may or may not vanish on the 
boundary. As such this equation cannot be affected by 
the addition of a boundary. To determine the boundary 
conditions we need to extend the support of the varia- 
tion to include the boundary. We have that the normal 
variation of the bulk energy can be written as 

8 x F b = J dA [FK$ + 2 T ab K ab <S> + T ab 8^K ab ] , (24) 

where r ah = dF/dj ab and T ab = dT/dK ab . The bound- 
ary term we wish to identify in 5j_F b originates in the 
8±K ab term in this expression. We recall that the ex- 
trinsic curvature transforms as follows under a normal 
deformation of £ (see e.g. M): 



8j_K ab 
We thus have that 



-V a V b <$> + K ac K c b <$>. 



S ± F b = / dA 



8$ 



'•( 



<$>V b T ab - T ab V b <5> 



(25) 



(26) 



where we have defined the Euler-Lagrange derivative 

£ = (- V a V 6 + K ac K c b )T ab + FK + 2 T ab K ab . (27) 

Thus, modulo the bulk shape equation, £ = 0, we get 
that the boundary contribution is 



S ± F b = j>^dsl a [§V b T ab - T ab V b <S>] 



(28) 



The terms proportional to V a $ and $ are not indepen- 
dent: the projection of V a $ along the edge is completely 
determined once $ is specified on C. To decompose S±F 
into two independent parts we proceed as follows: we first 
decompose V a $ into its normal and tangential parts with 
respect to C, 



V a $ = l a V_L$ + t„$ : 



(29) 



where we have defined = l a V a - We now perform 
an integration by parts on the $ term to obtain for the 
second term on the right hand side of Eq. (E8) 



ds l a T ab V b <S> 



ds 



l a l b T ab V - <$>—(l a T ab t b 



(30) 



where we have discarded a total derivative term with 
respect to arclength. In this way we succeed in isolating 
the independent normal variations at the boundary, the 
coefficients of $ and V^$. 

From Eqs. @, @, @, we obtain that the total 
boundary contribution of the normal variations is 



S± F b 



ds 



J b F ab V ± <S> 



tb 



aK 



}, (31) 



so that we can immediately read off the two boundary 
conditions which supplement Eq. (|2J 



laV b T ab + —(l a T ab t b 



+ aK\\ 
l a hF ab 



0, 
0. 



(32) 
(33) 



The first is of third order in derivatives of the embedding 
functions. This is consistent with the fact that the shape 
equation ( |l9| ) is of fourth order. Using the decomposition 
of the covariant derivative ( p9[ ) , it can be written in the 
alternative form 

l a l b V^T ab + 2^ s (l a T ab t b ^+K(lJ b -t a t b )T ab +aK\ ] = 0. 

(34) 

In the case of a membrane described by the Helfrich 
Hamiltonian (17) with an edge the third boundary con- 
dition ( |33"1) implies 



K = Kr 



(35) 



on the edge — the rigid membrane necessarily has a con- 
stant mean curvature at the edge equal to its spontaneous 
value. This is entirely independent of the tensions, /1 or 
a, or of the rigidity modulus. If Kq = 0, the membrane 
is minimal at its edge. 

The second boundary condition ( |3^ ) is of Robin type. 
For any T which is a function only of K, we have that 
T ab cx 7 ab , so that the middle term in (p2n vanishes, 



l a T ab t b — , 
and the boundary condition reduces to 



2a V 1 K 



aK\\ = . 



(36) 



(37) 



This equation determines the normal derivative of K in 
terms of the component of the extrinsic curvature tan- 
gent to the edge. It does not involve the surface tension 
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\x. We emphasize that its existence seems to have gone 
unnoticed. 

The first boundary condition, Eq. (p3|) , together with 
Eq. (p5[) , implies that on the edge 



<7K + fl = . 



(38) 



The geodesic curvature of a loaded boundary is com- 
pletely fixed by the ratio of the tensions in exactly the 
same way as in the previous section for soap bubbles, see 
Eq. ( |l2| ) . If fi = the edge is necessarily a geodesic of 
the bulk geometry. 

If the line tension on the boundary vanishes, a = 0, the 
consistency of Eq J3|) with ( p3| ) requires that fi = also. 
Furthermore Eq. (|37|) implies V ±K = on the boundary. 
But the unique solution satisfying the two boundary con- 
ditions K = Kq and W±K = is K = Kq everywhere. 
One way to see this is to construct the Gaussian nor- 
mal coordinates adapted to the edge, (I, s), where I is 
the length of the geodesic which intersects the edge nor- 
mally. With respect to this system of coordinates, the 
Laplacian assumes the form V 2 = df + ndi + d 2 in the 
neighborhood of the edge. Thus, modulo the boundary 
conditions, V 2 K = dfK on the edge. But Eq. (|l9|) im- 
plies that V 2 K = there so that dfK and all higher 
derivatives vanish. If K is analytic in I, then K = Kq. If 
p 0, there is no such constraint. The geometry is very 
severely constrained by the boundary conditions. 

Let us now examine an axially symmetric membrane 
with an axially symmetric edge. With respect to cylin- 
drical polar cooordinates {p, z, p} on i? 3 , the membrane 
is described by p = R(£), z = Z(t) where Z' 2 + R' 2 = 1. 



J 



£ is the arc length along a a curve with fixed ip, and the 
primes denote a derivative with respect to £. The intrin- 
sic geometry of E is described by the line element 



dT 2 



d£ 2 + R 2 (£)dp 2 



(39) 



We can write the extrinsic curvature in a form consistent 
with axial symmetry as 



K ab = £ a £ b K e + ( lab - U b )K f 



(40) 



where Ki and Kr are two spatial scalars which we iden- 
tify as the principal curvatures of the embedding of E in 
R 3 , £ a is the outward pointing unit normal to the circle of 
fixed £, £ a = (1, 0). The mean curvature is K = Ki+Kr. 
To evaluate the principal curvatures, it is convenient to 
define 8 as the angle which the tangent to a curve of 
fixed tp makes with the cylindrical radial direction: 



dZ 

— - = tan 6 . 
dR 

We then have Z 1 = sin0, and R' 
principal curvatures are 



(41) 



cosO, so that the 



K e = Q', 



K R = 



sin 9 
R 



(42) 



Axial symmetry implies that the fourth order shape equa- 
tion can be integrated to provide a third order equation 
for R as a function of I. It has been shown elsewhere 
(O, O], sec also 0) that this equation takes the form 



2a cos 9 9' 



sin 9 
R 



sin G 
R 



)(G' 



sinG N 
R ' 

I — 



) sinG + 2aK - 



sin 2 G 



R 



(fi + aKl) sinG = 0. 



(43) 



If the boundary C is also axially symmetric so that it 



coincides with a fixed value of I then l a 



Kr, K± = Ki and Kin — 0. It is simple to show that 
K = —R'/R. We thus have for the boundary conditions, 
Eqs.© and @, 

G' + sinG/i? = ^ , aR' = pR. (44) 

The remaining boundary condition, Eq. (|37|) , of third or- 
der in derivatives appears to present a problem: a third 
order ODE does not admit third order boundary con- 
ditions. The inconsistency, however, is only apparent: 
on the boundary, the shape equation Eq. ([43]) itself re- 
produces, modulo Eqs. (^4), the troublesome boundary 
condition Eq. (|37|) . Our analysis is thus completely con- 
sistent with the axially symmetric analysis of Q where 
the boundary conditions Eqs.(^) and (|3^) are derived. It 
is worth stressing, however, that potential pitfalls of us- 
ing the axially symmetric problem as a guide to the more 



general problem. The boundary condition (|3^) is a non- 
trivial constraint on the geometry which is not already 
encoded in the shape equation for non-axially symmetric 
configurations. 



IV. BALANCE OF FORCES 

In this section, we consider the balance of the forces 
operating at the edge. This provides the missing intuition 
on the physical origin of the boundary conditions we have 
derived in the previous section. 

Consider a point on the edge. In equilibrium, the ten- 
sion g must satisfy 



g 



(45) 

Here P is the membrane stress tensor so that f a l a is the 
'surface tension' acting on the edge due to unbalanced 
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stresses in the bulk at its boundary. In ||, it was shown 
that the bulk stress tensor for the model defined by the 
Helfrich energy ( |l7j ) can be expressed in the form 



f a = [2aK(K ab - y7 Q ' 



) - 2aK Q (K ab - K 1 ab ) 
2aV a Kn. (46) 



Thus its projection along the normal to the edge l a is 

f a l a = {2a(K-K )K ± -a(K-K ) 2 - f i}l 

+ 2a{K - K )K\ l± t-2aVxKn. (47) 



In addition, as we have seen in Sect. II, 



g = -<rt . 



(48) 



UsingEq. (|) for t, we read off the three components of 
Eq.@§, 

g k = 2a(K - K )K ± - a(K - K Q ) 2 - (i , (49) 
aK\\ = -2a\7 ± K, (50) 
= 2a(K-K )K l{± , (51) 



The condition ( |50| ) co- 
If K U + 0, 



respectively along 1, n, and t. 
incidcs with the boundary condition (37). 
Eq.(pll) implies that K = Kq. The remaining boundary 
condition ( pi| ) then coincides with a linear combination 
of the boundary conditions (pa) and ( 38 ) . In the axially 



symmetric geometry, however, Kin does vanish so that 
( j5l| ) does not imply K — Kq as it stands. One needs then 
to appeal to the integrated shape equation (f43|), which 
together with ( f49| ) and (|5^) reproduce K = K a . 

We thus have identified a very simple (if heavily dis- 
guised) physical interpretation of the boundary condi- 
tions. In particular, in this approach, the boundary con- 
dition K = Ko emerges as the vanishing of the stress 
induced by the bulk along the edge. Note that the vari- 
ational approach did not rely on the identification of 
projections. Indeed, the boundary condition correspond- 
ing to the projection along t was originally identified by 
demanding stationary energy for independent boundary 
variations of Vi$. 

We also note that the form of Eq.(|45|) implies the in- 
tegrability condition, 



dsf a l a = 



(52) 



on the edge. The existence of these three extremely non- 
trivial conditions is far from obvious in our previous ap- 
proach. 

One can say more. Take the equation, V a f a = 0, de- 
scribing the conservation of the stress tensor, dot it into 
X and integrate over the membrane surface. We get 



J dAV a (X-t a ) = JdAe a -t a . 



(53) 



Working on the right, we have 

JdAV a (X-t a ) = <j)dsX-lJ a 
dsX-g 



ds t ■ g 



(54) 



where we have used Eq.(|45|) on the second line, as well as 

Eq.(ff8|) on the last line. On the other hand, if we write 
p = jab Gh + then 



dAe n ■ i a 



dAf\. 



(55) 



We note that the bending energy J dA K 2 is a conformal 
invariant, and so does not contribute to the trace f a a . 
We have f a a = 2(aK a (K - K ) - fi), so that 

2fiA - 2aK J dA (K - K ) + oL = . (56) 

This condition is useful for identifying the sign associated 
with the mulupliers. For example, if Kq = 0, it is clear 
that fi is necessarily negative. 



V. GAUSSIAN RIGIDITY 

The geometrical scalars we can construct with dimen- 
sion [Length] ~ 2 are TZ, K 2 and K ab K ab . The Gauss- 
Codazzi equation ( p0| ) tells us that the three scalar in- 
variants TZ, K 2 and K a f,K ab are not independent. In 
addition, the Gauss-Bonnet functional 



1= / dATl 



(57) 



is a topological invariant if the membrane is closed. More 
generally for an open membrane, 



7= / dATl 



2 <j) ds k 



(58) 



is a topological invariant. A consequence is that if a 
Gaussian rigidity term is included in the energy a line 
rigidity § dsn is necessarily induced along its boundary. 

To obtain the variation of the Gaussian term, we need 
to know how the scalar curvature varies. Its tangential 
deformation is straightforward: to determine its normal 
deformation, we exploit Eq.(^0|) and with it the technol- 
ogy developed Sects. II and III. 

Consider now a Gaussian rigidity addition to the bulk 
energy, so that 



F = F b +f3 / dATl 



(59) 
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Whereas the bulk shape equation is unmodified, all three 
boundary conditions are changed: 



a(K - K ] 

2aV±K~ 20K\\_l + <tK 
a(K - K ) + 0K> 



(60) 
(61) 
(62) 



We note that the Gauss-Codazzi equation (gfj) allows us 
to express 1Z in terms of the projections of K a },_with 
respect to the edge, 1Z — 2(K^K± — K 



ll-L 



Eq.@ 



is 



quadratic in the extrinsic curvature. Eqs.(61) and (62j) 
by constrast are linear relationship between K± and Ku . 

Note that, unlike the case of the pure Helfrich model, 
the central term in Eq.(^Tj) does not vanish in general. 
However, it does vanish in an axially symmetric geometry 
(with axially symmetric edge), K\\j_ = 0. More generally, 
we have the integral statement 



j> ds [2a V^K + ctK{\ = 0. 



(63) 



In an axially symmetric geometry one can check that, 
modulo the lower order boundary conditions (|6^) and 
(|62l), Eq.(43) reproduces Eq.(|3l|) on the boundary. 

Let us consider now the balance of the forces in this 
case. The Gaussian term makes no contribution to f a 



Ij. Naively reinvoking Eq.(45) would appear to suggest 
that this term cannot modify the boundary conditions, in 
contradiction with what we have just derived. However, 
with a general function T and in particular for 01Z, the 
term d/ds(l a J 7ab £&)*£ appearing in its normal variation 
(see Eq.(|30[)) will be non- vanishing, and it is no longer 
appropriate to discard a total derivative as we did in 
deriving Eq. (|3C|). For consistency, we claim therefore 
that we need to modify Eq.([l5|) as follows: 



g - g - l a T ab t b n . 



(64) 



For Gauss-Bonnet, the second term reads — 2/3-fT|ij_n. 
This mysterious term is precisely the tension associated 
with the edge energy § c ds n. The projections of Eq.ji^) 
along 1, n and t respectively then read 

a [(K - K Q ) 2 - 2{K - K )K ± ] - 20Kf\ ± 



+ uk + /i = , (65) 

2 aVj_K - 2/3#||_l + aK\\ = , (66) 

2 K u [a{K - K ) + 0K\\\ =0, (67) 

where we have used the fact that, at the edge, 

n = .Kj|t + .K'||j L l. (68) 

As was the case in the previous section, these coin- 
cide with the boundary conditions (|6p|), (|6l|), (|62] ) when 

K\\± ¥= o. 

VI. CONCLUSIONS 



Whereas for a soap film, it is very simple to identify the 
forces operating on the edge, and so read off the boundary 
conditions on the bulk geometry, such an approach is less 
obvious for a membrane. However, we have demonstrated 
how simple geometrical and variational arguments may 
be exploited to derive the boundary conditions on the 
lipid membrane geometry. We have made no restrictive 
assumptions about the symmetry of the configuration. 
We then showed how these boundary conditions emerge 
from a balance of forces projected along a basis of vectors 
adapted to the edge. 

We plan at a future date to examine the axially sym- 
metric shape equation subject to the Gauss-Bonnet mod- 
ified boundary conditions. Matching observed configura- 
tions with solutions of the shape equation in this more 
general setting provides a method, at least in principle, 
for determining the values of the two bulk rigidity mod- 
uli. 
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